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ABSTRACT 
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The  development  of  a method  for  uniquely  characterizing  the 
many  different  symbols  which  are  used  to  represent  data  and  having  this 
characterization  general  enough  to  allow  for  minor  variations  in  the 
symbols  would  enable  a computer  to  obtain  information  directly  from 
written  text. 

The  characterization  consisted  of  obtaining  a Fourier  Transform 
representation  for  each  character  of  the  alphabet  and  a study  of  the 
harmonic  spectral  density  in  order  to  obtain  a pattern  by  which  the 
character  may  be  recognized.  The  pattern  by  necessity  must  be  .inique 
and  recognizable.  After  establishment  of  the  reference  alphabet,  each 
character  was  distorted  in  a like  manner  so  that  the  effects  of  the 
distortion  upon  all  the  harmonic  coefficients  could  be  seen.  The 
letters.  A,  E,  R,  and  S were  studied  in  greater  depth  and  the  results 
are  given 


The  characterizations' obtained  from  these  methods  show  a 
definite  uniqueness  of  representation  and  accordingly  a reference 

alphabet  was  characterized.  The  resolution  used  in  obtaining  the 

\ 

Fourier  Transform  representation  seemed  to  be  the  major  limiting 
factor  as  to  the  amount  of  distortion  acceptableT^The  most 
important  outgrowth  of  the  thesis  was  the  fact  that,  by  using  only 
the  resolution  which  was  used  in  the  characterization  of  the  listed — 
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Chapter  1 


INTRODUCTION 


Objectives 

The  ability  to  rapidly  perform  repetitive  calculations  and 
store  large  amounts  of  data  are  only  two  of  the  many  advantages  of  a 
computer.  Before  either  of  these  abilities  can  be  used,  the  machine 
must  first  receive  Its  commands  and  data  from  a set  of  rigidly  stan- 
dardized symbols  which  man  uses  as  his  communications  link  with  the 
machine.  These  symbols  must  follow  strictly  uniform  patterns  as  to 
height,  width,  and  relative  positioning  one  to  the  other.  Because  of 
the  need  for  this  rigid  standardization,  man  must  do  the  translating 
of  his  knowledge,  which  is  usually  in  the  ^orm  of  printed  text,  into 
symbols  recognizable  by  the  machine.  Although  man  may  use  the  same 
set  of  symbols  as  that  of  a machine,  he  restricts  his  symbols  to  very 
little  uniformity.  This  is  shown  best  by  noting  that  the  English 
language  has  but  one  alphabet  consisting  of  twenty  six  letters 
(symbols),  yet  there  are  almost  unlimited  variations  of  each  letter  used. 
However,  in  most  cases  no  matter  what  variations  are  used,  man  is  able 
to  distinguish  one  letter  from  another.  If  a computer  was  likewise 
able  to  distinguish  one  nonuniform  character  of  the  alphabet  from 
another,  then  such  a machine  could  obtain  information  directly  from 
text  written  in  almost  any  form  or  format. 

In  the  machine's  case,  as  in  the  human  case,  there  are 
limitations  in  the  amoemt  that  a symbol  can  be  varied  and  still  be 
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recognizable.  The  limitations  are  almost  totally  dependent  upon  how 
Intelligent  the  machine  is  or  put  another  way,  how  much  information 
about  a spacific  character  the  machine  has  access  too.  Since  a 
rigidly  standardized  alphabet  requires  the  minimum  amount  of 
information  necessary  to  distinguish  one  character  from  another,  this 
is  the  place  to  begin  the  characterization. 


Method  of  Procedure 

As  mentioned  before,  one  of  the  computer's  greatest 
advantages  was  its  ability  to  rapidly  perform  repetitive  calculations. 
Thus,  if  a mathematical  model  were  used  to  represent  each  character 
of  the  alphabet,  then  the  machine  could  be  used  to  obtain  this 
representation.  The  specific  math  model  used  was  important  only  to 
the  point  that  it  be  understandable,  its  computations  be  relatively 
short,  and  it  must  uniquely  represent  each  letter  or  symbol.  The 
uniqueness  of  representation  was  the  most  important  aspect  of  the 
mathematical  model. 

After  a suitable  method  of  modeling  was  established  for  a 
particular  reference  alphabet,  the  characterization  of  each  symbol 
was  tested  to  see  how  much  distortion  was  permissible  while  still 
having  the  representation  recognizable.  Because  of  the  many 
different  ways  a character  may  be  distorted,  an  extremely  comprehen- 
sive study  would  have  to  be  made  if  one  wishes  the  representations  to 
extend  any  further  than  the  printed  alphabet.  This  paper,  however, 
limited  itself  to  the  considerations  of  only  the  Future  (5)  type 
capital  letters.  This  set  of  characters  was  used  as  the  reference 
set  referred  to  in  this  paper  because  of  their  geometric  structure 
and  uniformity  of  stroke.  The  type  also  seemed  to  give  the  best 


accessible  example  of  a commonly  acceptable  alphabet  which  was 
readable  and  easily  reproduced  by  hand  or  machine. 
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Although  the  entire  reference  alphabet  is  distorted  in  only 
a few  ways,  four  specific  characters  are  chosen  and  studied  in 
greater  depth.  The  results  of  the  studies  are  given  in  the  appropri- 
ate section.  They  are  not  conclusive,  but  they  do  show  certain 
interesting  aspects  about  the  mathematical  modeling  and  characteriza- 
tion schemes  used.  For  example,  the  characteristic  patterns  obtained 
by  the  limited  resolution  are  unique,  easily  recognizable,  and 
seemingly  can  withstand  any  normal  distortion  which  occurs. 
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of  a particular  character  and  comparing  this  shape  to  his  reference 
set.  Admittedly,  this  is  an  oversimplification  of  man's  mental 
process,  but  for  the  purpose  intended  the  anology  will  suffice. 
Therefore,  the  problem  is  that  of  finding  a suitable  reference  set 
and  a method  of  characterizing  the  set  in  such  a way  that  each 
character  is  uniquely  represented. 

The  available  characterization  techniques  are  numerous, 
and  all  use  an  optical  scanner,  in  one  way  or  another,  to  obtain 
information  about  a character.  A few  of  the  more  successful  methods 
are  stroke  processing  (1),  the  use  of  a feature-template  preprocessor 
for  feature  extraction  (6),  classification  in  ad  hoc  fashion  based  on 
strokes,  concavities,  enclosures  and  blobs  (6),  and  the  Clemens' 
Technique  (6).  All  of  these  techniques  have  had  only  limited  success 
when  applied  to  nonuniform  characters;  therefore,  the  characterizations 
obtained  for  each  letter  were  not  uniquely  recognizable.  The  method 
used  in  this  paper  was  basically  a preprocessing  of  the  character  to 
obtain  certain  information,  formulating  a mathematical  model  using  the 
information  obtained  from  preprocessing,  and  characterization  based 
on  the  coefficients  of  a unique  model. 

Character  Preprocessing 

Each  charcter  referred  to  in  this  paper  was  preprocessed  by 
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hand.  However,  the  method  was  the  same  as  one  which  could  be 
employed  by  an  optical  scanner,  and  as  such,  the  preprocessing 
technique  will  be  explained  from  this  point  of  view. 
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The  purpose  of  preprocessing  was  to  change  a particular  charac- 
ter into  digital  data  and  have  this  data  in  a specific  format.  To  do 
this,  the  character  was  scanned  to  find  the  uppermost,  lowermost,  right 
most,  and  leftmost  points  in  order  to  set  the  character  boundaries. 
These  boundary  lines  consisted  only  of  horizontal  and  vertical  lines 
where  any  two  adjacent  lines  form  a right  angle.  This  had  the  effect 
of  normalizing  the  letter  to  a rectangular  area.  The  rectangle  was 
then  subdivided  into  64  smaller  but  equal  areas  as  shown  in  Figure  1A 
and  these  areas  were  called  cells. 

To  obtain  digital  information,  the  cells  were  assigned  values 
of  zero  or  one  depending  on  the  contents  of  each  cell.  If  a cell 
contained  any  part  of  the  character,  it  was  assigned  a value  of  one; 
otherwise,  the  cell  had  a value  of  zero.  The  process  is  shown  in 
Figure  2A,  and  the  resultant  was  an  eight  by  eight  matrix  of  zeroes 
and  ones. 

Next,  the  matrix  was  read,  and  the  points  were  plotted  as  dis- 
crete samples  with  uniform  spacing  on  an  x-y  graph  (rectangular 
coordinate  graph).  Certain  letter  characteristics  could  be  observed 
by  the  particular  choice  of  a starting  point  and  direction  of  reading. 
After  deliberation  over  the  different  symmetrical  characteristics  of 
each  letter  of  the  alphabet,  a method  was  chosen  so  that  the  symmetry 
of  a reference  character  about  a vertical  line  drawn  through  the  center 
of  the  matrix  could  be  observed.  The  method  of  reading  consisted  of 
starting  at  ceil  (1,4)  which  was  the  first  point,  and  reading  toward 
the  left  until  cell  (1,1)  was  reached.  Then  the  same  procedure  was 
repeated  from  cell  (2,4)  to  cell  (2,1)  and  repeatedly  until  cell  (8,1) 


was  reached.  These  points  were  then  plotted  parallel  to  the  negative  x 
axis  of  the  x-y  graph.  The  next  points  to  be  read  were  cells  (1,5)  to 
(1,8).  Then,  repeating  this  left  to  right  movement,  cells  (2,5)  to  (2,8) 
were  read.  This  was  done  until  cell  (8,8)  had  been  reached.  These 
points  were  plotted  parallel  to  the  positive  x axis  on  the  x-y  graph 
with  the  first  point  being  cell  (1,5). 

The  result  of  this  method  of  reading  was  that  it  showed  a char- 
acter which  was  symmetric  about  a vertical  line  drawn  through  the  center 
of  the  matrix  as  having  half-wave  symmetry.  However,  the  characteriza- 
tion of  a letter  with  such  symmetry  gave  a large  d.  c.  component.  To 
avoid  this,  the  half-wave  symmetry  was  changed  to  negative  halfnove 
symmetry  by  changing  the  sign  on  the  points  located  in  the  left  half 
plane  of  the  x-y  graph.  In  order  to  better  visualize  the  shape  of  the 
graph,  each  sample  point  was  connected  by  a straight  line  to  the  next 
one.  Then  the  resultant  waveform  was  shifted  32  sample  points  to  the 
right  so  as  to  locate  the  entire  wave  in  the  right  half  plane  as  shown 
in  Figure  3A. 

The  final  preprocessing  step  was  obtaining  data  from  the  x-y 
graph  which  will  be  used  in  calculating  the  mathematical  model  of  a 
letter.  This  was  accomplished  by  digitally  sampling  the  waveform  at 
128  equally  spaced  intervals.  The  sign  of  each  data  point  must  be 
preserved. 

Mathematical  Modeling 

The  character,  at  this  point,  was  represented  by  a particular 
waveform  on  an  x-y  graph.  This  indicated  that  the  character's  x-y  plot 


Pattern  Recognition 

Because  of  the  unique  Fourier  representation  for  each  reference 
letter,  characterization  based  upon  the  coefficients  of  this  representa 
tion  was  the  final  step.  All  characterizations  obtained  in  this  paper 
were  achieved  by  manual  means.  Computer  oriented  inethods  are  possible 
and  would  result  in  faster  and  more  exact  characterizations. 

The  manual  method  was  as  follows.  After  obtaining  the 
coefficients  of  the  Fourier  representation  for  a particular  letter,  the 
magnitudes  of  the  coefficients  at  each  harmonic  were  calculated  as 

(1*)  IRMS*2{k)l  = J (S(k))^  + (T(k))^  ‘ . 


These  magnitudes  were  plotted  versus  the  harmonic  number  as  a bar 
graph.  A bar  graph  of  reference  letter  A is  shown  in  Figure  4A. 

This  same  process  was  used  for  each  reference  letter  and  also  for  each 
distorted  letter.  Next,  all  of  the  plots  for  a particular  letter  were 
compared  to  find  consistent  similarities.  These  similarities  can  be 
mathematically  defined  for  single  harmonics  as 
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Figure  4A  Sample  Bar  Graph  ( Magnitude  Spectrum  ) 


13 


where  IDCI  is  the  coefficient  inagnitude  of  the  distorted  character  at  a 
particular  harmonic,  IRCI  is  the  coefficient  magnitude  of  the  reference 
letter  at  a specific  harmonic,  and  IRCFI  is  the  coefficient  magnitude 
of  the  reference  letter's  fundamental  harmonic.  In  the  case  of  group 
harmonics,  the  defining  equation  is 


(6) 


IHMDGl  - IHMRGI  ^ 
I RFC  I 


^ Difference  for  Group  Harmonics 


where  IHMDCI  is  harmonic  mean  of  a particular  group  in  the  distorted 
symbol's  magnitude  spectrum  and  IHMRGI  is  the  harmonic  mean  of  a 
particular  group  in  the  reference  letter's  spectrum.  The  harmonic 
mean  can  be  calculated  as 


(7) 


IHHI  - 


N,  + 1) 


where  and  are  respectively  the  left  and  rightmost  harmonics  of 
the  group  and  Xj  is  magnitude  of  the  ith  harmonic.  These  similarities 
formed  the  characterization  of  the  letter.  A precaution  that  no  two 
reference  characters  have  the  same  representation  should  be  taken. 

This  was  especially  true  for  the  manual  method.  As  described,  the 
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manual  method  had  many  undesirable  aspects.  Hovw'ever , it  served 
adequately  for  the  purposes  of  this  paper. 


1i» 


A superior  method  wouid  have  been  a computer  oriented  approach. 
Such  a method  of  pattern  recognition  would  be  a necessity  if  the 
reading  of  nonuniform  characters  were  to  be  handled  by  a computer.  The 
advantages  of  using  a computer  would  be  in  several  ways.  First,  a 
matrix  of  higher  resolution  (having  more  than  Sk  cells)  could  be  used 
which  would  mean  that  the  fine  details  of  a particular  character  could 
be  represented.  Also,  the  characterization  of  a particular  reference 
letter  could  be  carried  out  to  far  greater  extremes  which  would  decrease 
the  possibility  of  having  the  same  representation  for  two  reference 
letters. 

There  are  many  techniques  already  developed  that  will  suffice 
if  the  plot  (coefficient  magnitude  versus  harmonic)  was  considered  as 
a picture.  If  such  a consideration  was  made  and  since  the  plot  will 
always  be  of  the  same  nature,  certain  techniques  could  be  applied  such 
as  a least  mean-square  error  fit  (2),  derivative  data  feature 
extraction  (3),  use  of  a contextual  analysis  and  topological  coding 
program  (4),  or  an  extension  of  the  manual  method  used  in  this  paper. 
Possibly  a combination  of  several  methods  would  be  the  most  advantageous 
approach. 


Chapter  3 

DETAILED  EXAMINATION  OF  FOUR  CHARACTERS 

The  hand-printed  characters  were  formed  without  the  use  of 
a straight  or  guidelines. 

Character  A 

The  character  A was  subjected  to  four  types  of  distortion. 

The  reference  letter  A is  shown  in  Figure  5A  while  the  four  distorted 
A characters  are  shown  in  Figures  6A,  7A,  8A,  and  9A.  Comparison 
of  the  distorted  characters  to  the  reference  symbol  in  a numerical 
sequence  showed  a steady  increase  in  the  deviation  of  the  characters 
from  the  reference.  A similar  sequential  look  at  Figures  5B  through 
9B  gave  more  insight  into  the  effects  of  various  distortions.  Several 
interesting  points  should  be  noted  at  this  time.  First,  the  width  of 
the  stroke  used  to  form  a letter  seemed  to  have  very  little  effect 
upon  the  matrix  except  in  the  cases  of  extremely  broad  or  narrow  widths 
where  these  might  cause  parts  of  the  letter  to  become  obscured. 
Secondly,  visual  inspection  of  the  matrices  showed  that  although 
character  symmetry  was  lost  the  letter  was  still  easily  recognizable 
in  roost  cases.  Also,  at  this  point  it  should  be  noted  that  in  the 
norroal ization  of  script  character  A,  the  head  and  foot  finals  (5) 
have  been  eliminated.  By  thinking  of  these  strokes  as  only  letter 
connections  in  script-written  text,  the  elimination  of  them  was 
Justified  because  they  added  little  or  nothing  to  the  character 
recognition  process.  Upon  close  inspection  of  the  script  A matrix 
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In  Figure  9B,  one  would  also  notice  that  the  letter  which  this  matrix 
defined  was  not  recognizable  within  the  matrix  as  the  other  A characters 
were. 

A sequential  comparison  of  each  character's  matrix  with  the 
resulting  waveform  was  next.  Figures  5C  through  9C  are  the  waveshapes. 
From  this  comparison,  the  symmetry  of  the  A reference  character  about 
a vertical  line  drawn  through  the  center  of  its  matrix  could  be  observed. 
Another  observation  made  was  the  fact  that  as  a letter  was  shifted 
more  to  one  side,  this  same  side  of  the  matrix  had  an  increase  in  the 
ratio  of  valued  cells  to  zero  cells.  This  increase  could  be  seen  in  the 
waveshape  as  a shifting  of  area.  Next,  by  looking  only  at  the  waveform 
of  each  of  the  characters,  a distinct  similarity  between  most  of  them 
could  be  seen. 

However,  the  heart  of  this  analysis  lay  in  the  magnitude 
spectrum  of  each  of  the  characters.  Figures  50  througi  90  show  the 
bar  graphs.  As  stated  earlier,  the  characterization  of  a reference 
letter  was  contained  in  the  similarities  of  these  magnitude  spectra. 

There  were  no  indepth  computer  comparisons  made  of  the  bar  graphs 
because  of  the  limited  distortions  and  with  only  visual ly-obvious 
similarities  being  noted.  These  similarities  included  the  strict 
magnitude  at  a specific  harmonic,  the  magnitude  relative  to  nearest 
neighbors  at  a specific  harmonic,  combinations  of  strict  and  relative 
magnitude  changes  in  a group,  and  general  shape  of  the  bar  graph. 

The  characterization  of  reference  letter  A was  made  without 


regard  to  script  character  A.  The  reason  for  disregarding  the  A script 


symbol  could  best  be  seen  in  a comparison  of  its  matrix  and  wave- 
shape to  the  matrix  and  waveshape  of  the  reference  character.  Although 
both  characters  were  known  to  be  the  letter  A,  the  matrix  for  the  script 
A did  not  have  enough  resolution  to  allow  the  symbol  to  be  recognizable 
in  it.  This  being  the  case,  the  Fourier  representation  of  the  matrix 
could  not  be  used  in  forming  the  characterization  for  the  reference  A. 

An  increase  in  resolution  could  be  accomplished  only  if 
several  facts  were  kept  in  mind.  The  first  one  was  that  the  number 
of  positive  and  negative  lobes  in  the  waveform  was  directly  related 
to  the  number  of  rows  in  the  matrix,  and  increasing  the  number  of  rows 
would  negate  any  useful  characterization  information  already  obtained. 

So  the  increase  in  resolution  should  be  accomplished  by  increasing 
the  number  of  columns  which  only  increases  the  number  of  data  points 
and  the  emphasis  on  small  details  in  the  character.  However,  because 
of  the  severe  distortion  in  the  script  A character,  such  an  increase  in 
the  number  of  columns  still  could  not  catch  enough  detail  to  allow 
character  recognition  based  on  the  reference  A.  In  order  to  obtain 
character  recognition  of  such  a letter,  a minimum  matrix  of  32  by  32 
would  have  to  be  used  in  characterizing  the  reference  set. 

For  reference  symbol  A,  the  following  characterization  was 
made.  Harmonic  numbers  15,  27,  30,  3^,  and  49  were  chosen  as  the 
similar  single  harmonics.  By  referring  to  Table  1A,  the  largest  percent 
difference  could  be  seen  to  be  almost  17  percent  which  was  for 
harmonic  15  of  the  slanted  A and  the  hand-printed  A.  Harmonic  numbers 
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Figure  5D  Magnitude  Spectrum  for  Reference  Character  A 
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Figure  6A  Light  Character  A 
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Figure  6B  Matrix  for  Light  Character  A 
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Ftgure  6D  Magnitude  Spectrum  for  Light  Character  A 
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Figure  8D  Magnitude  Spectrum  for  Hand-Written  Character  A 
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Figure  9D  Magnitude  Spectrum  for  Script  Character  A 


Table  1A 

Percent  Difference  between  the  Characterizing 
Single  Harmonics  of  the  Distorted 
and  Reference  A Letters 


IRC  I Percent  Difference 


15  0.04  0.36  0.36  0.08  0.46  10.0  16.7  16.7  63.3 

27  0.02  0.02  0.04  0.04  X 2.5  2.5  5.8  5.8 

30  0.04  0.06  0.04  0.06  0.02  3.3  6.7  3-3  6.7 


34  X 0.08  0.06  0.02  0.02  2.5 


6.7  0.0 


49  0.04  0.04  0.04  X 0.06  3.3  3.3  3.3  9.2 


Column  1 refers  to  the  Light  character 
Column  2 refers  to  the  Slanted  Character 
Column  3 refers  to  the  Hand-Printed  Character 
Column  4 refers  to  the  Script  Character 
X implies  a Bar  Graph  value  of  zero  but  for  calculation 
purposes  a value  of  .005  is  assumed 


Table  IB 
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Percent  Difference  between  the  Characterizing 
Group  Harmonics  of  the  Distorted 
and  Reference  A Letters 
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1 through  3,  18  through  21,  4l  through  kj , and  5^  through  63  were 
chosen  as  the  similar  group  harmonics.  By  referring  to  Table  IB,  the 
largest  percent  difference  could  be  seen  to  be  almost  17  percent  which 


was  for  the  group  1 through  3 of  the  hand-printed  A.  Since  no  indepth 
analysis  was  obtained,  the  visual  appearance  of  the  bar  graphs  was  very 
important.  By  t'  criteria  the  magnitude  spectra  of  letters  showed  a 
large  amount  of  resemblance. 

Character  E 

The  character  E was  subjected  to  four  types  of  distortion.  The 
reference  character  E is  shown  In  Figure  10A,  and  the  four  distorted  E 
characters  are  shown  in  figures  11A,  12A,  13A,  and  14A.  A comparison  in 
numerical  sequence  of  the  reference  character  to  the  distorted  letters 
showed  a steady  increase  in  the  deviation  of  the  distorted  characters 
from  the  reference.  A similar  sequential  look  at  Figures  10B  through 
14B  allowed  a better  understanding  of  the  effects  of  various  distortions. 
At  this  time  several  points  should  be  noted.  First,  the  stroke  width 
used  to  form  a letter  had  little  or  no  effect  on  the  matrix  except  in 
the  extreme  cases  noted  for  character  A.  Secondly,  visual  inspection 
of  the  matrices  showed  that  although  the  physical  geometry  of  a letter 
was  altered,  the  character  was  still  recognizable.  It  should  also 
be  noted  that  there  was  no  need  for  letter  truncation  as  experienced 
with  character  A.  Because  of  this  and  the  fact  that  the  script 
character  E was  recognizable  in  its  defining  matrix,  the  script 
character  was  used  in  the  character  recognition  of  the  reference  set. 
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A sequential  comparison  of  each  character's  matrix  with  the 
resulting  waveform  was  next.  Figures  IOC  through  14C  are  the  wave- 
shapes. From  this  comparison,  a distinct  feature  of  printed  characters 
as  opposed  to  script  characters  was  observed.  The  printed  character 
showed  its  basic  use  of  straight  vertical  and  horizontal  lines  as  single 
large  areas  of  the  x-y  graph  where  as  the  curved  formations  of  the 
script  E showed  up  as  smaller,  more  numerous  areas,  provided  the 
resolution  of  the  matrix  was  adequate. 

However,  the  bar  graph  of  each  character  more  fully  described 
the  outcome.  Figures  10D  through  14D  show  the  magnitude  spectra.  As 
stated  earlier,  the  characterization  of  a reference  character  was 
determined  by  the  similarities  of  its  bar  graph  and  those  of  the  dis- 
torted characters.  The  same  criteria  was  used  for  this  ietter  as  was 
used  for  character  A.  ^ 

I 

1 

The  following  characterization  for  reference  letter  E was  made.  \ 

j 

Harmonic  numbers  2,  5,  11,  26,  28,  3k,  37  and  39  were  chosen  as  the  j 

similar  single  harmonics.  By  referring  to  Table  2A,  the  largest 
percent  difference  could  be  seen  to  be  almost  32  percent  which  was  for 
the  fifth  harmonic  of  the  script  E.  This  large  percent  difference 
seemed  to  be  caused  by  the  head  and  foot  finals  oi  the  script  E. 

Harmonic  numbers  7 through  9,  17  and  18,  3^  through  40,  and  44  through 
63  were  chosen  as  the  similar  group  harmonics.  By  referring  to  Tabie 
2B,  the  largest  percent  difference  could  be  seen  to  be  almost  28  per- 
cent which  was  for  the  group  7 through  9 of  the  light  E.  An 
examination  of  the  magnitude  spectra  of  the  light  and  reference  E 
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Figure  10D  Magnitude  Spectrum  for  Reference  Character  E 
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Figure  12B  Matrix  for  Slanted  Character  E 
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Figure  13D  Magnitude  Spectrum  for  Hand-Written  Character  E 
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Table  2A 

Percent  Difference  between  the  Characterizing 
Single  Harmonics  of  the  Distorted 
and  Reference  E Letters 
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Column  1 refers  to  the  Light  Character 
Column  2 refers  to  the  Slanted  Character 
Column  3 refers  to  the  Hand-Printed  Character 
Column  4 refers  to  the  Script  Character 
X implies  a Bar  Graph  value  of  zero  but  for  calculation 
purposes  a value  of  .005  is  assumed 


kk 


symbols  showed  that  the  targe  difference  was  due  exclusively  to  the 
eighth  harmonic  of  the  light  E.  However,  because  there  was  such  a great 
visual  resemblance  between  the  bar  graphs  of  the  letters,  the  character- 
ization was  not  changed. 

Character  R 
' ' 

Character  R was  subjected  to  four  types  of  distortion.  The  j 

reference  character  R is  shown  in  Figures  15A  while  the  four  distorted  || 

’] 

R letters  are  shown  in  Figures  16A,  17A,  18A,  and  I9A.  By  comparing  | 

i 

the  distorted  letters  to  the  reference  symbol  in  a numerical  sequence, 
an  increase  in  deviation  of  each  character  from  the  reference  was  seen. 

A similar  comparison  of  Figures  15B  through  19B  gave  a better  under- 
standing of  the  effects  caused  by  various  distortions.  | 

Several  interesting  points  could  be  noted  at  this  time.  First, 
as  with  the  two  previous  characters,  varying  the  width  of  the  stroke  || 

used  to  form  a character  seemed  to  have  very  little  effect  on  the 
defining  matrix  unless  the  width  was  carried  to  an  extreme  as  in  Figure 
17B.  Secondly,  as  the  letter  was  slanted  to  one  side,  there  was  a 
resultant  shift  in  area  of  the  waveshape.  At  this  point  it  should  be 
noted  that  in  the  normalization  of  script  symbol  R,  the  head  and  foot 
finals  (5)  have  been  eliminated  as  was  the  case  for  script  letter  A. 

The  same  justifications  that  were  used  for  script  A were  used  here 
for  script  R character  truncation.  Also,  upon  inspection  of  the  matrix 
in  Figure  19B,  one  would  note  that  the  letter  which  this  matrix 
defined  was  not  recognizable  within  the  matrix. 

A sequential  comparison  of  each  character's  matrix  with  the 
resulting  waveshape  was  next.  Figures  15C  through  19C  are  the 
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Figure  16A  Light  Character  R 
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Figure  16B  Matrix  for  Light  Character  R 
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Figure  16D  Magnitude  Spectrum  for  Light  Character  R 


Figure  I7A  Slanted  Character  R 
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Figure  170  Magnitude  Spectrum  for  Slanted  Character  R 


Figure  18D  Magnitude  Spectrum  for  Hand-Written  Character  R 
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Figure  19D  Magnitude  Spectrum  for  Script  Character  R 
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Table  3A 

Percent  Difference  between  the  Characterizing 
Single  Harmonics  of  the  Distorted 
and  Reference  R Letters 


IRC  I Percent  Difference 


0.20  0.06  0.10  0.20  0.12  8.9  6.7  2.2 

0.08  0.08  0.12  0.06  0.10  2.2  2.2  2.2 


30  0.02  0.02  0.02  0.04  0.02 
40  0.02  0.02  X 0.04  0.02 
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Column  1 refers  to  the  Light  Character 
Column  2 refers  to  the  Slanted  Character 
Column  3 refers  to  the  Hand-Printed  Character 
Column  4 refers  to  the  Script  Character 
X implies  a Bar  Graph  value  of  zero  but  for  calculation 
purposes  a value  of  .005  is  assumed 
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waveforms.  In  this  comparison  the  distortion  due  to  a large  increase 
in  the  width  of  the  stroke  was  very  apparent  in  Figure  17C. 

However,  the  slant  R character  was  still  clearly  visible  in  its 
defining  matrix.  But  in  the  case  of  the  script  R,  a comparison  showed 
that  the  matrix  was  unable  to  define  the  letter,  and  thus,  its 
magnitude  spectrum  was  not  used  in  the  characterization. 

Figures  150  through  190  show  the  bar  graphs  and  foIloi.ing  the 
same  guidelines  as  used  previously,  the  characterization  of  reference 
letter  R was  determined.  Harmonic  numbers  2,  S,  30,  and  kO  were 
chosen  as  the  similar  single  harmonics.  The  largest  percent  difference 
as  seen  from  Table  3A  was  almost  9 percent  which  was  for  the  second 
harmonic  of  the  light  R character.  Harmonic  numbers  3 through  6,  16 
through  19*  2k  through  27,  k3  through  46,  and  48  through  63  were  chosen 
as  the  similar  group  harmonics.  The  largest  percent  difference  was 
almost  7 percent  as  shown  in  Table  3B  for  the  16  through  19  group  of 
harmonics  of  the  slanted  R letter.  An  interesting  point  should  be 
noted  at  this  time  about  the  percent  differences  of  the  script  R 
character,  in  most  all  cases  its  difference  was  larger  respectively 
than  for  any  other  letter.  This  showed  that  the  magnitude  spectrum 
of  the  script  R was  not  one  which  could  be  recognized  as  representing  an 
R. 

Character  S 

Reference  letter  S and  the  four  distorted  S characters  are 
shown  in  Figures  20A  through  24A.  A comparison  in  numerical  sequence 
of  the  distorted  letters  to  the  reference  symbol  showed  a steady 
increase  In  the  deviation  of  each  distorted  character  from  the 
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reference.  A similar  review  of  Figures  20B  through  24B  showed  the 
effects  of  each  distortion. 

Here  again  the  stroke  width  used  seemed  to  have  little  effect 
on  the  matrix  except  in  extreme  cases  such  as  Figure  22B.  Although  the 
slant  S character  was  clearly  distinguishable  as  the  equivalent  of  the 
reference  S symbol,  the  slant  S letter  was  barely  definable  in  its 
matrix.  Although  this  was  the  case,  the  slant  S's  magnitude  spectrum 
was  still  used  in  the  characterization  of  the  reference  S letter. 

Visual  inspection  of  the  other  matrices  showed  that  with  limited 
geometrical  letter  distortion  the  other  S characters  could  still  be 
recognized  in  their  defining  matrix.  When  comparing  Figures  2kA 
and  24b,  it  should  be  noted  that  no  symbol  truncation  was  necessary 
for  the  script  character  S.  Because  of  the  fact  that  the  script  S 
was  recognizable  in  its  defining  matrix,  the  symbol  was  used  in  the 
characterization  of  the  reference  letter. 

A comparison  of  each  character's  matrix  with  the  resulting 
waveform  was  next.  The  waveshapes  are  pictured  in  Figures  20C 
through  24C.  This  comparison  showed  the  nature  of  curved  symbols  as 
opposed  to  the  straight  line  composed  letters  viewed  until  now.  The 
semicircular  composition  of  the  S was  not  recognizable  in  the  wave- 
shapes so  even  if  an  S were  composed  of  straight  horizontal  and 
vertical  lines,  it  would  still  appear  the  same.  This  means  that  a 
symbol  did  not  have  to  be  a smooth  evenly  contoured  stroke;  instead, 
it  could  consist  of  short  continuous  straight  lines  and  still  be 
read 1 1 y recogn i zed . 


5S 

For  the  characterization  of  the  S reference  ietter,  Figures 
200  through  240,  which  show  the  bar  graphs,  was  used.  Also,  the  same 
criteria  which  was  used  for  the  previous  characters  was  used  for  this 
symbol . 

For  reference  symbol  S,  the  foiiowing  characterization  was 
made.  Harmonic  numbers  2,  9,  21,  24,  32,  and  40  were  chosen  as  the 
similar  singie  harmonics.  The  largest  percent  difference  as  shown  in 
Table  4A  was  slightly  more  than  10  percent.  This  occurred  for  the 
ninth  harmonic  of  the  script  S character.  Harmonic  numbers  3 through 
5,  11  and  12,  29  through  31,  and  44  through  63  were  chosen  as  the 
similar  group  harmonics.  The  largest  percent  difference  as  shown  in 
Table  4B  was  almost  26  percent  which  was  for  the  3 through  5 group  of 
harmonics  for  the  script  S letter.  It  was  interesting  to  note  at  this 
time  that  the  greatest  percent  differences  did  not  occur  in  the 
slanted  S letter,  although  this  particular  symbol  was  barely  recogniz- 
able in  its  defining  matrix.  Also,  the  overall  visual  appearance  of 
the  distorted  and  reference  letter's  bar  graphs  have  a close  resem- 
blance. 
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Figure  20D  Magnitude  Spectrum  for  Reference  Character  S 
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Figure  210  Magnitude  Spectrum  for  Light  Character  S 
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Figure  22D  Magnitude  Spectrum  for  Slanted  Character  S 
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Figure  23D  Magnitude  Spectrum  for  Hand-Written  Character  S 


Figure  2kD  Magnitude  Spectrum  for  Script  Character  S 
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Table  ^tA 

Percent  Difference  between  the  Characterizing 
Single  Harmonics  of  the  Distorted 
and  Reference  S Letters 


IRCI  Percent  Difference 


Column  1 refers  to  the  Light  Character 
Column  2 refers  to  the  Slanted  Character 
Column  3 refers  to  the  Hand-Printed  Character 
Column  4 refers  to  the  Script 

X implies  a Bar  Graph  value  of  zero  but  for  calculation 
purposes  a vaule  of  .005  is  assumed 


I 
I 
I 

I Table  kB 

Percent  Difference  between  the  Characterizing 
f Group  Harmonics  of  the  Distorted  I 

and  Reference  S Letters  ; 


Har. 
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Chapter  k 


SUMMARY 


Analysis  of  Results 

When  analyzing  the  results,  the  main  objective  of  this  paper 
should  not  be  forgotten.  Although  the  characterization  of  a reference 
alphabet  was  frequently  mentioned  and  used,  this  was  done  oniy  as  a 
means  to  an  end  and  not  as  a final  result.  This  means  that  the 
establishment  of  characterized  reference  letters  were  used  as  a 
measure  of  the  process,  and  the  actual  end  result  was  whether  or  not 
the  techniques  used  would  give  unique  representations. 

Because  of  the  Fourier  Transform  properties,  there  should  be 
no  doubt  that  the  representation  obtained  for  a particular  normalized 
matrix  was  unique.  However,  because  of  the  limited  resolution  used  in 
a matrix,  it  would  be  possible  to  have  two  different  characters  with 
equivalent  matrices  which  would  result  in  identical  magnitude  spectra. 
So  the  first  weak  point  in  the  process  as  outlined  by  this  thesis  was 
the  degree  of  resolution  necessary  to  avoid  duplicate  matrices  for 
different  letters.  Since  the  ideal  case  of  infinite  resolution  was  not 
practical,  the  decision  of  how  much  resolution  to  use  was  based  on  the 
limitations  placed  on  the  process.  Of  these  limitations  the  two  most 
important  were  the  time  required  to  character  recognize  a symbol  and 
the  maximum  distortion  of  any  letter. 
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These  two  requirements  were  not  independent  of  each  other;  in 
fact  one  was  directly  related  to  the  other.  This  was  easily  understood 
because  as  a symbol's  distortion  increases  so  does  the  resolution 
necessary  to  define  the  character.  The  results  were  more  data  points 
to  be  transformed  which  took  longer  computer  time. 

There  was  only  one  other  point  in  the  recognition  scheme  that 
could  invalidate  this  system.  This  weak  point  was  in  the  method  used 
to  characterize  a reference  letter.  The  characterization  of  reference 
symbol  was  based  on  the  similarities  between  its  magnitude  spectrum 
and  those  of  the  same  letter  but  with  distortions  introduced.  For  the 
purposes  of  this  thesis,  visual  means  for  defining  the  similarities 
sufficed,  but  for  a working  system  better  recognition  methods  would  be 
needed.  Several  possibi 1 ? ties  were  discussed  previously;  however,  along 
with  these  methods  the  maximum  allowable  distortion  must  again  be 
considered  before  setting  error  criteria. 


a 


Any  further  investigation  into  the  procedure  outlined  in  this 
thesis  should  be  directed  toward  the  recognition  process.  Although 
two  weak  points  were  discussed,  the  problem  of  matrix  resolution  was 
only  a difficulty  for  the  particular  considerations  of  this  paper. 

The  resolutions  which  are  obtainable  by  optical  scanning  are  more  than 
adequate  for  a working  system.  So  the  need  for  an  automated 
recognition  process  is  paramount,  and  one  of  the  major  considerations 
of  such  a process  would  be  control  of  the  error  parameters.  This 
control  would  be  ncessary  in  order  to  obtain  a valid  and  unique 
representation  for  each  member  of  the  reference  set. 


Although  the  technique,  as  described  by 
for  character  recognition  of  distorted  symbols. 


this  thesis,  did  allow 
it  did  not  completely 


solve  the  other  problems  encountered  by  an  actual  working  system. 


1 

1 

1 

I 


I 

4 


A: 

Appendix  A 

°0 

Reference  Set  I 

This  appendix  contains  those  reference  characters  not 
previously  discussed.  Their  magnitude  spectra,  which  resulted  from  | 

the  computer  program  of  Appendix  B,  are  also  Included. 
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Figure  26A  Reference  Character  C 
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Figure  27A  Reference  Character  D 
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Figure  28A  Reference  Character  F 
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Figure  30A  Reference  Character  H 


Figure  37A  Reference  Character  0 


Figure  38A  Reference  Character  P 


Fjgure  39A  Reference  Character  Q 


Figure  hOA  Reference  Character  T 


Figure  41A  Reference  Character  U 


Figure  42A  Reference  Character  V 


Figure  43A  Reference  Character  W 


Figure  45A  Reference  Character  Y 


Figure  46A  Reference  Character  Z 
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Figure  28B  Hagnitude  Spectrum  for  Reference  Character  F 
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Figure  31B  Magnitude  Spectrum  for  Reference  Character  I 
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Figure  32B  Magnitude  Spectrum  for  Reference  Cherecter  J 
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Figure  33B  Magnitude  Spectrum  for  Reference  Character  K 
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Figure  34B  Magnitude  Spectrum  for  Reference  Character  L 
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Figure  35B  Magnitude  Spectrum  for  Reference  Character  H 
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Figure  36B  Magnitude  Spectrum  for  Reference  Character  N 
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Figure  38B  Magnitude  Spectrum  for  Reference  Character  P 
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Figure  IiOB  Magnitude  Spectrum  for  Reference  Character  T 
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Figure  41B  Hagnitude  Spectrum  for  Reference  Character  U 
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Figure  426  Magnitude  Spectrum  for  Reference  Character  V 
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Figure  436  Magnitude  Spectrum  for  Reference  Character  W 
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Figure  kSB  Magnitude  Spectrum  for  Reference  Character  Y 


Figure  Magnitude  Spectrum  for  Reference  Character  Z 


Appendix  B 

Computer  Program 

The  following  Is  an  explanation  In  the  use  of  the  computer 
program  that  was  used  for  the  purpose  of  obtaining  data  necessary  for 
this  thesis.  A listing  of  the  program  wr I tten  in  Fortran  IV 
Immediately  follows  this  explanation.  Since  the  program  is  divided 
into  the  main  program  and  two  subroutines,  each  part  Is  discussed 
individual ly. 

Main  Program 

The  main  program  basically  consisted  of  the  calling  sequences 
for  the  two  subroutines,  a manipuitation  of  data  into  the  correct  form 
for  use  by  the  subroutines,  and  the  Input/output  operations.  Because 
the  input  has  the  possibility  of  being  complex  numbers,  the  read  format 
statemnt  must  take  this  into  account.  Therefore,  the  Input  should  be 
broken  into  its  real  and  imaginary  parts.  The  real  parts  of  the  data 
were  placed  first  on  the  cards  according  to  the  format.  Then  the 
imaginary  parts  of  the  data  followed  according  to  the  same  format. 

This  method  of  reading  placed  the  real  and  imaginary  parts  of  each  data 
point  In  adjacent  memory  locations  which  was  the  necessary  arrange- 
ment  for  the  F.  F.  T.  subroutine. 

Manipulation  of  data  by  the  main  program  changed  the  output 
of  the  F.  F.  T.  subroutine  into  complex  numbers  which  were  the 
coefficients  of  the  sine  and  cosine  terms  in  the  Fourier  expansion. 
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Then,  these  coefficients  were  used  to  calculate  a number  which  was 
proportional  to  the  total  root-mean-square  value  at  each  harmonic. 

This  number  was  to  be  used  by  the  plot  routine. 

F.  F.  T.  Subroutine 

The  basis  for  the  algorithm  of  the  F.  F.  T.  subroutine  was  the 
composite  nature  of  certain  integer  numbers  and  the  cyclic  properties 
of  wJk 

where  W is  represented  as 

(8)  W « exp  ( — and  i - 


Hathematically,  the  F.  F.  T.  subroutine  was  derived  from  the  Discrete 
Fourier  Transform  (D.F.T.)  whose  forward  transform  can  be  represented 
as 


(9)  .j  - x(k)  . 1^'  x(k)  wJ" 

k-0  k^D 

where  J * 0,1,2, N-1 

This  equation  Is  already  In  an  algorithmic  form  which  can  be  readily 
programmed  as  a single  program  on  a digital  computer. 

By  placing  a limitation  on  N (the  number  of  complex  data  points) 
which  required  It  to  be  a power  of  two,  the  cyclic  nature  of  W-^^,  as 
mentioned  previously,  became  apparent  In  that  its  only  possible  values 
were  ♦•1,  -1,  +1,  and  -I.  In  general  repeated  with  every  four 


I 

» consecutive  Integer  numbers.  Thus,  by  prearranging  the  data  in  a 

specific  order,  this  repetition  could  decrease  the  number  of  calcula- 


tions necessary  for  a transform.  Specifically,  the  number  of  complex 
multiplications  and  additions  necessary  in  the  F.  F.  T.  subroutine 
contained  in  this  thesis  is 


(10) 


N log2  N (7) 


where  N is  the  number  of  complex  data  points  and  a power  of  two.  The 
savings  provided  by  this  subroutine  were  very  substantial  for  large 

i 

values  of  N.  For  example,  the  OFT  would  require  over  a million  complex 
multiplications  to  transform  1,024  data  points,  but  the  F.  F.  T. 
subroutine  required  about  5>000  complex  multiplications.  This  was  a 
reduction  of  about  200  to  1. 

Besides  performing  the  actual  transformation,  this  subroutine 
checked  N to  assure  that  it  was  a power  of  two,  and  also  preshuffled 
the  input  data  into  the  correct  sequence  as  required  by  the  F.  F.  T. 
algorithm.  The  preshuffling  was  strictly  based  on  N and  the  four- 
integer-repeti tion  cycle  of 

Plot  Subroutine 

The  plotting  subroutine  takes  information  already  contained  In 


memory  and  plots  these  values  against  the  harmonic  number.  The  scale 
factors  are  also  set  in  this  subroutine. 
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INTEGER  X1,Y1,W1,W2,W3 
COMPLEX  RK,SK,TK,VI 
N-128 
W1-N/2 

DIMENSION  A(2,128)  ,X(61t)  ,ARMS2(64) 

READ(5,1) ((A(X1,Y1) ,Y1-1 ,128) ,X1-1 ,2) 

1 FORMAT(8F5.1) 

CALL  FFT(A,N,-1) 

CON- 1. /FLOAT (N) 

WRITE(6.100) 

100  FORMATdHI  ,19X,1HJ,AX,8HRE(A(J))  ,l4X,8HIM(A(J))//) 

DO  2 1-1 ,N 
K-l-1 

AR-A(1,I)*C0N 

AI-A(2,I)*C0N 

2 WRITE(6,110)K,AR,AI 

110  F0RMAT(18X,I3,2(1*X,F8.5)) 

WRITE(6,112) 

112  F0RMAT(1H1,11X,1HK,13X,AHR(K) ,25X,^HS(K) ,25X,AHT(K)  , 
*15X,8HRMS(K)*2,5X,//) 

DO  95  1=1 ,W1 
Vl-(0,1) 

K-l-1 

L-N-K+1 

IF(K.EQ.0)AA-(A(1 ,l))*CON 

IF(K.EQ.0)BB-(A(2.I)*C0N 

IF(K.GT.0)AA=(A(1,I)+A(1 ,L))*CON 

IF(K.GT.0)BB-(A(2,I)+A(2,L))*C0N 

IF(K.GT.0)DD=(A(1,I)-A(1,L))*C0N 

IF(K.GT.0)EE-(A(2,I)-A(2,L))*C0N 

IF(K.E(l.O)RK-CMPLX(AA.BB) 

IF(K.EI1.0)SK-(0,0) 

IF(K.EQ.0)TK-(0.0) 

IF(K.GT.0)RK-(0.0) 

IF(K.GT.O)SK-CMPLX(AA.BB) 
IF(K.GT.O)TK-(CMPLX(DD,EE))*VI 
ARMS2 1-CABS (CSQRT( (SK**2)+(TK**2) ) ) 
ARMS2(I)-ARMS21’A=1000.0 
95  WRITE(6, 111)K,RK,SK,TK,ARMS21 

111  F0RMAT(10X,13,7(4X,FT0.6)) 

DO  94  1-1 ,W1 

94  X(l)-FLOAT(l) 

WRITE(6,113) 

1 13  FORMAT('l') 

CALL  PL0T(X,ARMS2) 

86  STOP 
END 
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SUBROUTINE  FFT(D,N,NSIGN) 

COMPLEX  T,D(N) ,W,WARG 
TEST  THAT  INPUT  IS  POWER  OF  TWO 
1-1 

1 l-l+l 
IF(I-N)1,3,2 

2 WRITE(6,100) 

100  FORMAT ( IX, 40HDATA  BLOCK  NOT  POWER  OF  TWO.  CAN  NOT  RUN) 
RETURN 

START  OF  DATA  SHUFFLE 

3 NHF-N/2 
JO-O 

LAST-N-2 
DO  8 I -2, LAST 
M-NHF 

k IF(JO.LT.M)GO  TO  5 
JO-JO-M 
M-M/2 
GO  TO  A 
5 JO-JO+M 
J-JO+1 

IF(I .GE.J)GO  TO  8 
T=D(J) 

D(J)-D(I) 

D(I)=T 

8 CONTINUE 
BEGIN  FFT 
INC-1 
MNODES-2 

9 IF(MNODES.GT.N}  RETURN 
DO  10  M-1,INC,1 

WARG-CMPLX (0 . , 6 . 283 1 853* FLOAT (NS  I GN* (M- 1 ) ) /FLOAT (MNODES ) ) 
W-CEXP(WARG) 

DO  10  l-M,N, MNODES 
J-I+INC 
T-0(J)*W 
D(J)-D(I)-T 
10  D(l)-D{l)+T 
INC-MNODES 

MNODES-MNODES+MNODES 

GO  TO  9 

END 
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SUBROUTINE  PLOT(X  ,ARMS2) 

XMAX-70.0 
XHIN-0.0 
YMAX»1 200.0 
YMIN-0.0 
DATA  IB.IA/' 

DIMENSION  IP1(70)  ,1X1(8)  ,IP(120,70)  ,X{64)  ,ARMS2(6l*) 
1-120 

YMAX 1 - ( YMAX- YM I N/FLOAT ( I ) 

DO  92  L-1,120 
DO  2 K»1,70 

2 IP(L,K)=IB 

92  IP(L,K)-IB 

5 DO  3 J-1,6I» 

IY1-(ARMS2(J)-YMIN)/YMAX1 
IF(IY1.NF.I)G0  TO  7 
IF(IYl.EQ.l) lY-lYl 
IF(lY1.E(i.l)G0  TO  8 

7 IY-(ARMS2(J)-YMIN)/YMAX1+1.Q 
IF(IY.NE.I)G0  TO  3 

8 IX-(X(J)-XMIN-1 
IX3-IX 

DO  93  13-1, lY 

93  IP(I3,IX3)-IA 

3 CONTINUE 
11-1 

DO  91  12-1,70 
1X2-12 

91  IP1(I2)-IP(I1 ,1X2) 

IF((l/10)*10.E(l.l)G0  TO  4 
WRITE(6,20) IP1 

20  FORMATdIX,'*' ,70A1) 

GO  TO  6 

k YP-(YMAX1/1000)*I+YMIN 
WRITE(6,21)YP,IP1 

21  F0RMAT(1X,F8.1,2X,'+',70A1 

6 1-1-2 

IF(I .GT.0)G0  TO  5 

YP-YMIN 

WRITE(6,22)YP 

22  F0RMAT(1X,F8.1,2X,'+' ,7 ( '*********+')) 
XUNIT-((XMAX-XMIN)/70.)*10. 

DO  9 1-1,8 

9 IX1(I)-XUNIT*(I-1)+XMIN 
WRITE(6, 12) (1X1(1) ,1-1,8) 

12  F0RMAT(2X,8(7X,I3)) 

WRITE(6,120) 

120  F0RMAT('0' ,10X,'PL0T  OF  RMS*2  VS  HARMONIC  NUMBER') 
WRITE(6,11A) 
nit  FORMAT('I') 

RETURN 

END 
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